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We present a model for unicellular algal growth as motivated by several experiments 
implicating the importance of calcium ions and "loosening" enzymes in morphogenesis. A 
growing cell at rest in a diffusive calcium solution is viewed as an elastic shell on short 
timescales. For a given turgor pressure, we calculate the stressed shapes of the wall elements 
whose elastic properties are determined by Young's modulus and the thickness of the wall. 
The local enzyme concentration then determines the rate at which the unstressed shape 
of a wall element relaxes toward its stressed shape. The local wall thickness is calculated 
from the calcium-mediated addition of material and thinning due to elongation. We use 
this model to calculate growth rates for small perturbations to a circular cell. We find an 
instability related to modulations of the wall thickness, leading to growth rates which peak 
at a finite wave number. 
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In recent years, there has been increasing interest in models of unicellular growth. Specific examples include 
dendritic branching in neuronal growth and the broadening and branching of lobes in unicellular algal 
growth (H ||]. The algal models have ranged from purely geometrical models (analogous to the geometrical 
models of dendritic growth in solidification to one which combines the geometric formalism with a 
diffusive "morphogen" field. At the same time, a wealth of experimental information regarding unicellular 
morphogenesis has been provided by studies of the larger species of the alga Micrasterias. Morphogenesis in 
Micrasterias proceeds by a very well-ordered sequence of tip-splitting and lobe-broadening, culminating in 
elegant fan-shapes (Fig. I) which are reminiscent of the patterns seen in unstable diffusion-controlled growth 
j7j . This paper is devoted to the construction of a new model for this process, a model built on the dynamics 
of the cell wall. 

Although some of the details remain elusive, a general picture of cell growth in algae has emerged from 
experiment. In order for the cell wall to elongate without thinning indefinitely, vesicles containing the 
CLi appropriate building materials are synthesized within the cytoplasm and then travel toward the cell periphery 

where the fusion of vesicles with the cell membrane is thought to be mediated by calcium ions J§|]. In 
particular, growing tips of Micrasterias exhibit high concentrations of membrane-associated calcium [^|,^), 
enabling them to fuse vesicles in much larger quantities than other areas of the cell wall JlQ] . Meanwhile, 



there is evidence that calcium concentration does not vary significantly within the body of the cell itself |L6| , 
perhaps implying that the relevant diffusive processes occur outside the cell body even as they modulate the 
concentration at the cell wall. 

Experiments that halt growth by reducing turgor pressure demonstrate that elongation is (at least in part) 
a response to the stresses in the cell wall [jlCj . In addition, it has long been believed that a "loosening factor" 
must be present to allow the fibers making up the wall to slip past each other during elongation (e.g. the 
protein "expansin" pj]|). Since most elongation occurs at the cell tips p 12 1, we might surmise that high 



concentrations of calcium also imply high concentrations of the "loosening" factor. These observations have 
led us to a simple model for cell wall dynamics in algae. 

Consider a cell growing in a solution of calcium ions. We treat the cell wall as an elastic shell on short time 
scales, whose slow plastic deformation is then catalyzed by the concentration of a "loosening" factor. We 
further assume that the fusion of vesicles with the cell membrane (and subsequent discharge of the vesicle 
contents into the cell wall) is also a long time scale process. 

We begin by discretizing the stressed shape into N rods with lengths U and orientations . Likewise, the 
unstressed shape will consist of N rods with lengths l[ and orientations 9[ (and our convention hereafter will 
be to associate all primed quantities with the unstressed shape). The stressed shape minimizes the energy 
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FIG. 1. Micrasterias denticulata, diameter ~ 200/im. Picture provided by U. Meindl. 
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where the angle made by the intersection of two successive rods gives rise to the curvatures Ki ~ 
(Qi+i — 0%) j\ (h+i + h)- The sums are, respectively, the energy due to the pressure (p § xdy, where p < 
for an outward pressure), the usual strain energy arising from Hooke's Law (where at = Ehi/l[ and E is 
Young's modulus), and the pure bending energy of a rod (where = EI^ and the moment I = J x 2 dx in- 
tegrated over the rod's width for an isotropic material Jl^,0 ) ■ Note that if the cell has no known rotational 
symmetry, we must also constrain the cell shape to be closed. 

Once the stressed shape is determined, we can quasistatically calculate the calcium concentration at 
the cell boundary. As suggested by experiments, we assume that diffusion occurs outside the cell, with 
the concentration C{R oa ) = Coo. We therefore solve V 2 u = in the region outside the cell, where u = 
(C — Coo) /Coo is the normalized concentration and u{Roo] — 0. At the cell boundary itself, we specify the 
flux through the membrane with du/dn — j(u)/D, where n is the outward normal to the cell surface, j(u) 
is the inward current, and D is the diffusion constant. 

We then allow the unstressed shapes to relax toward the stressed configuration via 

dY 

-±=T{u){k-l[) (2) 



d0 l 

^A=T(u)(9 i -9' i } (3) 

where T(u) is the characteristic relaxation rate as a function of the "loosening" factor assuming, for simplicity, 
that it is the equivalent of the calcium concentration. Finally, if we let x ( u ) be the rate at which mass 
(equivalently, the area hil'j, if the density is uniform) is added to a given rod through the fusion of vesicles, 
then 

dh % = x (u) hi dl'j , . 

dt I'. l'.dt' U 

which completes the specification of the model. 
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We can now use this model to calculate the growth rates of perturbations to the cell shape. As the 
physics of the problem is invariant under reparametrizations, we can always choose to parametrize the 
unstressed membrane by dividing it into N rods of equal length l\ — a . Then 9[ = n/2 + fa + b[S sin mfa 



- mb' 1 R'^~ 1 8 cos m<f) where 
. We expect the stressed 
[Note: Hereafter, we will 



with fa = 2ir (i — 1) /N, specifies a perturbed circle with a curvature k' w R' ~ 
R' = Na' /2ir. We specify the perturbed wall thickness as hi = ho + hiScosmq 
coordinates to have the form k — ciq + a±6 cos mfa and 6i = n/2 + fa + b\S sinmfa 
refer to the pieces of a function F expanded in 5 by F = Fq + Fi8 cos mfa] 

In the continuum limit, we expand the total energy as U « Uq + S 2 U2, noting that the O(S) energy vanishes 
upon integration with respect to <fi fl7|| . We solve the O (<5°) problem by minimizing Uq with respect to ao, 
yielding a quartic equation for ao = 2ttRq/N. 
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EhoRQ + EIqRq — EIqR'q = 



(5) 



While we could obtain the exact solution for ao by taking the appropriate root, it is more enlightening to 
expand the solution to O (p), which gives i?o ~ R'o —pR'o I (EhoR' +EIo). we can solve for the perturbations 
by minimizing U2 with respect to a\ and b\ (a long but straightforward calculation), yielding an answer in 
terms of ao |20|] . 
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where B = ERl(I R' ~ h §^R ) - EI (Rq - R' )(Ro ~ R'o§r)> J o = H h o), and dl/dh is the derivative 
evaluated at h = ho- Observe that if the bending energy is made to vanish for a cell with uniform thickness 
(i.e. I = dl/dho = hi = 0) we find that the stressed shape is a perfect circle (ai =61 =0). This agrees 
with the well-known "membrane" result that the tension T = \p\/n, implying that curvature variations in a 
pliable membrane require external support Jig) ]. 

We must now calculate the concentration on the stressed shape, which by direct integration of dxi — lj cos 6i 
and dyi = h sin 6*; is seen to have the radial perturbation Ri = Ro (mbi — ai/ao) / (to 2 — 1) for to > 1. The 
unperturbed solution is easily found to be uq (r) — RoD~ 1 jo ln(r/i?oo) where jo = j(uo(Ro))- Meanwhile, 
the perturbed solution must have the form u(r,fa) « uo (r) + cr ~ m S cos m<j). Applying the flux boundary 
condition, we find the concentration at the cell wall, u(R) = u + uiS cosmfa to be p2[ 



u(R) » u (Ro) + Jo( - m — ^RiS cosmfa (8) 



where dj/duo (and like expressions will) refer(s) to the derivative evaluated at u — uq(Ro)- 

The stressed and unstressed shapes, along with the concentration just calculated, allow us to calculate 
l[ (t + dt) and 9[ (t + dt) through equations (||) and (||). Reparameterizing the solution at t + dt so that we 
again have equal length rods yields 

dt ma Q ouo to 

where To = T (uq (Ro))- The rate of wall-thickening can be calculated directly from (0) as 



^ = - + [u! + ho^do) + T Q (Wo + ^0^) j S COSTO0 (10) 

dt a' n { \ a' Q duo ou Q J \ a' Q J J 
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where xo = X ( u o (Ro)) an d £o = ( a o — a d)/ a d- Since the average wall thickness is observed to be relatively 
constant during the growth of the cell 0, we must choose xo = a oha^o^o- We have tested these analytic 
results against simulations of the instantaneous rates and found them to be in agreement, with the error 
converging as 1/N 2 . 

Though these equations are somewhat complicated, we can still use them to examine qualitative behavior 
in several simple situations. For definiteness, we assume that both the relaxation rate and the vesicle fusion 
rate increase with higher concentrations (i.e. dT/duo, dx/duo > 0). The bending moment is taken as / oc h n 
where n > 1 and /i<i? |l4|]. Also, in accordance with experimental observations B, we assume j (u) causes 
an inward current at the cell tips (e.g. j = u — c u with c u a constant). For later reference, note that the 
stressed curvature is n « Rq 1 + Rq 1 {mb\ — ai/ao) Scosmcp. 

First, consider a perfectly circular cell that develops develops a slight thickening of the cell wall at <j> = 
(i.e. b[ — and hi > 0). For small pressures we find that (to lowest order) b\ = — £o^i(-fo — h n -^)/mI > 
and a\ = —hi(ao — a' )/ho < 0. This implies that tf> = is a "tip" (ki > 0), which is also a minimum of 
strain (a\ < 0). This agrees with the experimental observation that wall stresses are minimized at cell tips 
p3| . Since u± > at tips, Tohoai/ao is the only stabilizing term in (|Io|). Therefore, if ho\ai\/ao is not too 
large, the thickness variations will grow, providing positive feedback to the instability which sprouted this 
tip. 

Another interesting case is that of a perturbed cell shape which has a constant thickness (i.e. b[ > 
and hi = 0). To lowest order in pressure, we find that {b\ — b[) ~ P^Rq /Eh^ 1 (m 2 — l) < 0, a\ ~ 
— (ao — a' ) mb'^hQ /Rq < 0, and m ~ vnb\. Again, we find that the minimum strain occurs at the "tip". 
Looking at equation (^), we see that — T$ai/ao is destabilizing, while toFo(oi — 6^) is stabilizing. But since 
h <c i?o, we expect the net effect to stabilize b[. Apparently, without modulations in the cell wall thickness, 
tips will be smoothed out. ft is worth noting that this effect may explain the "lobe- broadening" observed in 
later stages of tip growth, though only full numerical simulations would demonstrate this. 

A trivial yet interesting implication of equation ( |ltj| ) is the ability to reproduce the observed patterns of 
deposition of wall material when the turgor pressure is reduced jL0|. Setting p = implies that £ = a\ = 0, 
eliminating all stabilizing influences in the thickness and implies dh\/dt = Ui Jj^-/a . That is, the thickness 
variations will simply follow the variations in concentration, allowing large amounts of material to collect at 
tips. 

Finally, we can use equations @ and (|l0|) to obtain a pair of equations of the form dbi/dt = Cn&j + C12/11 
and dhi/dt = ci\b' x + C2ih\. By assuming the coefficients vary slowly and letting b[, hi ~ exp(Xt), we can 
obtain an expression for the quasi-static growth rates. While it is unprofitable to write the expression here, 
we nonetheless plot the result for a typical set of parameters in Figure 2. We find that small-scale distur- 
bances are damped out, with the growth rate peaking at a finite wave number determined by the physical 
parameters of the model. In addition, we note that the m = 1 perturbation is indeed a "zero-mode." 

Spectrum of Growth Rates 

0.00006 - 




FIG. 2. Spectrum of eigenvalues in the quasistatic approximation. Growth rates (A) are plotted against wave 
number (m) for p = -0.001, E = 20.0, h = 0.4, R' = 10.0, R x = 10 6 , D = 1, j = (u + 0.15), V = u + 0.15, 
X = 0.12 (u + 0.15). 
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It is interesting to see that this relatively simple model can exhibit such rich behavior and reproduce 
several experimentally observed effects. This model has explicity assumed that reshaping of the cell wall is 
a relaxational process, wherein turgor pressure deforms the wall while enzymes allow the wall elements to 
slowly assume these stressed forms as their permanent forms. Our model provides a way to understand these 
long-standing biological assumptions analytically. It also forms the physical basis for numerical simulations 
that may provide further insight into the mechanisms of morphogenesis at later stages of growth. 

This work has been supported in part by NSF Grant No. DMR94-15460, and by the San Diego chapter 
of the ARCS Foundation. 
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